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Abstract 

This work is devoted to non-linear stochastic Schrodinger equations with mul- 
tiplicative fractional noise, where the stochastic integral is defined following the 
Riemann-Stieljes approach of Zahle. Under the assumptions that the initial con- 
dition is in the Sobolev space H q (W l ) for a dimension n less than three and q an 
integer greater or equal to zero, that the noise is a Q— fractional Brownian motion 
with Hurst index H S (^,1) and spatial regularity H q+4 (M. n ), as well as appro- 
priate hypotheses on the non-linearity, we obtain the local existence of a unique 
pathwise solution in C°(0, T, H g (M. n )) n C 0,7 (0, T, H q ~ 2 (R n )), for any 7 6 [0,H). 
Contrary to the parabolic case, standard fixed point techniques based on the mild 
formulation of the SPDE cannot be directly used because of the weak smoothing 
in time properties of the Schrodinger semigroup. We follow here a different route 
and our proof relies on a change of phase that removes the noise and leads to a 
Schrodinger equation with a magnetic potential that is not differentiable in time. 



1 Introduction 



This work is concerned with the existence theory for the stochastic Schrodinger equation 
with fractional multiplicative random noise of the form 

f = iAVdt -i^dBf 1 -ig(V)dt, t > 0, ieK", n<3 
I ¥(f = 0, .) = *„, 

where Bf = B H (t, x) is an infinite dimensional fractional Brownian motion in time and 
smooth in the space variables. The sense of the stochastic integral will be precised later 
on and the term g(^) is non- linear. We limit ourselves to n < 3 for physical considera- 
tions, but the theory should hold for arbitrary n with adjustments of some hypotheses. 
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Our interest for such a problem is motivated by the study of the propagation of paraxial 
waves in random media that are both strongly oscillating and slowly decorrelating in 
the variable associated to the distance of propagation. Such media are encountered for 
instance in turbulent atmosphere or in the earth's crust [HI 123]. More precisely, it is 
well-known that the wave equation reduces in the paraxial approximation [27] to the 
Schrodinger equation on the enveloppe function \1>, which reads in three dimensions 

% d z m = -A ± V + V(z, x)i&, z G M + , x e 1R 2 , 

where z is the direction of propagation of the collimated beam, x = (£1,0:2) is the 
transverse plane, Aj_ = d 2 2 + d 2 2 , and V is a random potential accounting for the 

x l x 2 

fluctuations of the refraction index. If V is stationary and its correlation function 
R(z, x) := E{V(z + u, x + y)V(u, y)} has the property that 

R(z,x) ~ z- a R (x), 

z—¥oo 

where Rq is a smooth function and < a < 1, then the process V presents long-range 
correlations in the z variable since R is not integrable. Rescaling V as V — > e~ zV(z/e, x) 
and invoking the non-central limit theorem, see e.g. [28], one may expect formally when 
e — > that 

-4-W-, x)^ £ (z, x) { ^ Wz, x)dB H (z, x), 

£2 e 

where B H is a Gaussian process with correlation function 

E{B B (z, x + y)B H (z', x)} = _ (z™ + {z'f H - \z - zf H )Ro(x), (2) 

with H = 1 — f G (5, 1). Proving this fact is an open problem, while the short-range 
case (when R is integrable) was addressed in [HEGI, and the limiting wave function \1/ is 
shown to be a solution to the Ito-Schrodinger equation. Our starting point here is (CQ), 
where we added a non-linear term g(^) to account for possible non-linearities arising 
for instance in non-linear optics. 

Let us be more precise now about the nature of the stochastic integral in ([T|). Since 
(JT]) is obtained after formal asymptotic limit of a L 2 norm preserving Schrodinger equa- 
tion, one may legitimately expect the limiting equation to also preserve the L 2 norm. 
The appropriate stochastic integral should therefore be of Stratonovich type, which in 
the context of fractional Brownian motions are encountered in the literature as pathwise 
integrals of various types, e.g. symmetric, forward, or backward, [2j [32] . Since in our 
case of interest the Hurst index H is greater than |, all integrals are equivalent and can 
be seen as Riemann-Stieljes integrals of appropriate functions, see [21 [32] and section [2] 
for more details. Such an integral is well-defined for instance if both integrands are of 
Holder regularity with respective indices (3 and 7 such that (3 + 7 > 1 [32] • 

In the context of SPDEs, the infinite dimensional character of the Gaussian process is 
usually addressed within two frameworks, whether for standard or fractional Brownian 
motions: the Q— (fractional) Brownian type, or the cylindrical type, see [1]. The first 
class is more restrictive and requires the correlation operator Q in the space variables 
to be a positive trace class operator (or even more for fractional Brownian motions, see 
[17]); in the second class, it is only supposed that Q is a positive self-adjoint operator on 
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some Hilbert space with appropriate Hilbert-Schmidt embeddings. As was done in [T7] 
for parabolic equations with multiplicative fractional noise, we will assume our noise is of 
Q— fractional type, which yields direct pathwise (almost sure) estimates on B H in some 
functional spaces. The cylindrical case is more difficult and our approach does not seem 
to generalize to it. The Q— fractional case actually excludes stationary in x correlation 
functions of the form (J2J) since they lead to a cylindrical type noise, which is a drawback 
of our assumptions. This latter case, even in the more favorable situation of parabolic 
equations, seems to still be open. Standard Brownian motions are more amenable to 
cylindrical noises since the Ito isometry holds. In the case of fractional type integrals, 
the "Ito isometry" involves the Malliavin derivative of the process, which is difficult to 
handle in the context of SPDEs with multiplicative noise. Hence, an existence theory 
for the Schrodinger equation in some average sense seems more involved to achieve, and 
we thus focus on a pathwise theory which requires the Q— Brownian assumption in our 
setting. 

Stochastic ODEs with fractional Brownian motion were investigated in great gener- 
ality in [2TJ . Stochastic PDEs with fractional multiplicative noise are somewhat difficult 
to study and to the best of our knowledge, the most advanced results in the field are 
that of Maslowski and al [17] . Duncan et al [10] or Grecksch et al [H]. The reference 
[TU] involves finite dimensional fractional noises, which is a limitation. Several other 
works deal with additive noise, which is a much more tractable situation as stochastic 
integrals are seen as Wiener integrals [TTJ [29j 03] an d cylindrical type noises are allowed. 
References [T7J [TUJ [H] consider variations of parabolic equations of the form 



where A is the generator of an analytic semigroup S, and the equation can be comple- 
mented with non-linear terms and a time-dependency in A in [TU]. The noise B H in 
these references is a Q— fractional Brownian with possibly additional assumptions. The 
difficulty is naturally to make sense of the term udBf 1 and to show that u is Holder 
in time. In that respect, the analyticity hypothesis is crucial: indeed, the standard 
technique to analyze fl3j) is to use mild solutions of the form 



and for the integral to exist, one needs the term S(t — s)u(s) to be roughly of Holder 
regularity in time with index greater than 1 — H . This means that both u and the 
semigroup S need such a regularity. While the term u(s) can be treated in the fixed 
point procedure, the semigroup S(t — s) has to be sufficiently smooth in time, which 
holds for analytic semigroups, but not in the case of Cq unitary groups generated by iA 
in the Schrodinger equation. In the latter situation, one can "trade" some regularity in 
time for S with some spatial regularity on u, but this procedure does not seem to be 
exploitable in a fixed point procedure. Another possibility could be to take advantage of 
the regularizing properties of the Schrodinger semigroup that provide a gain of almost 
half a derivative in space, and therefore to almost a quarter of a derivative in time [3J. 
It looked to us rather delicate to follow such an approach since the smoothing effects 
hold for particular topologies involving spatial weights which looked fairly intricate to 



du = Audi + udB 



t ■ 



(3) 
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handle in our problem, even by using the classical exchange regularity /decay for the 
Schrodinger equation. The strictly linear case (i.e. when g = in (PQ)) can likely 
be treated by somewhat brute force with iterated Wiener integrals and the Hu-Meyer 
formula, but this approach does not carry on to the non-linear setting. 

We propose in this work a different route than the mild formulation and a quite 
simple remedy based on two direct observations: (i) the usual change of variables for- 
mula holds for the pathwise stochastic integral and (ii) using it along with a change 
of phase removes the noise and leads to a Schrodinger equation with magnetic vector 
potential A(t,x) = —VB(t,x). Forgetting for the moment the non-linear term g{^), 
introducing ip(t,x) := e lB C ,:c )\I/(t, x) the filtered wavefunction, and supposing without 
lack of generality that B H (0,x) = 0, this yields the system 

p(* = 0, •) = *o, U 

which is a standard Schrodinger equation with a time- dependent Hamiltonian. There 
is a vast literature on the subject, see [211 13 US CSl HS1 UHl EH M 120] for a non- 
exhaustive list. One of the most classical assumptions on A for the existence of an 
evolution operator generated by the Hamiltonian A b h is that A is a C 1 function in 
time with values in H 1 (M n ). This is of course not verified for the fractional Brownian 
motion. The price to pay for that is to require additional spatial regularity, and one 
possibility (likely not optimal) is to suppose that B H has values in if 4 (R n ). Assuming 
such a strong regularity is naturally a drawback in this approach. 

Regarding the treatment of the non-linearity, we suppose that it is invariant by a 
change of phase, that is e lBt g(^f) = g((p), which is verified by power non-linearities of 
the form g(^) = or by g(^f) = V[^f]^ where V is the Poisson potential. Con- 

trary to the case of non-linear Ito-Schrodinger equations where various tools such as 
Strichartz estimates or Morawetz estimates have been successfully used to investigate 
focusing/defocusing phenomena in random and deterministic settings [31 [8j [7J [6] , there 
are very few available techniques to study (jlj) with a potential vector A not smooth in 
time and augmented with the term g. There are Strichartz estimates in the context of 
magnetic Schrodinger equations, but some require A to be C 1 in time [31] . and some 
others avoid such an hypothesis but assume instead that A is small in some sense [26J, 
which has no reason to hold here. As a result, we are lead to make rather crude as- 
sumptions on g in order to obtain a local existence result. Moreover, the analysis of 
non-linear Schrodinger equations generally relies in a crucial manner on energy methods. 
In our problem of interest, we are only able to obtain energy conservation for smooth 
solutions, which turns out to be of no use when trying to obtain a global-in-time result 
and limits us to local results, unless the non-linearity is globally Lipschitz in the appro- 
priate topology. This is due to the fractional noise that does not allow us to obtain H l 
estimates for \1/ via the energy relation, as we explain further in remark [2J 

The main result of the paper is therefore a local existence result of pathwise solutions 
to ([T]) with a smooth Q— fractional noise B B and appropriate assumptions on the non- 
linearity g. The article is structured as follows: in section [21 we recall basic results on 
fractional stochastic integration, and present our main result in section [3j section H] is 
devoted to the magnetic Schrodinger equation (jlj), while section concerns the proof of 
our main theorem. 
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2 Preliminaries 



Notation. We denote by H k (R n ) and W k ' q (R n ), 1 < n < 3, the standard Sobolev 
spaces with the convention that H°(M. n ) := L 2 (M n ). For a Banach space V, T > 0, 
and < a < 1, W a ,i(0, T, V) denote the space or mesurable functions / : [0,T] — > V 
equipped with the norm 

The space C°' a (0, T, V) denotes the classical Holder space of functions with values in V. 
When V = C or M, we will simply use the notations W a> i(0,T), C°' a (0,T) and || ■ || a i. 
Notice that for any e > 0, C°- Q+£ (0, T, V") C W a>1 (0, T, 1/).' For two Banach spaces U and 
V, C(U, V) denotes the space of bounded operators from U to V, with the convention 
C(U) = £(U, U). The L 2 inner product is denoted by (/, g) = f Rn fgdx where / is the 
complex conjugate of /. 

Fractional Brownian motion. For some positive time T, we denote by f3 = 
{(3 H (t), t G [0, T]} a standard fractional Brownian motion (fBm) over a probability 
space (f2,J-*, P) with Hurst index H G We will denote by L 2 (Q) the space of 

square integrable random variables for the measure P and will often omit the depen- 
dence of /3 H on u G Q for simplicity. The process /3 H is a centered Gaussian process 
with covariance 

E{/?f /3f } = l -(t 2H + s 2H -\t- s\ 2H ). 

Since E{(/3^ — P^) 2 } = \t — s| 2K , (3 H admits a Holder continuous version with index 
strictly less than H. In order to definite the infinite dimensional noise B H (t, x), consider 
a sequence of independent fBm (^) ngN . Let Q be a positive trace class operator on 
L 2 (IR n ) and denote by (fi n , e re ) n6 N its spectral elements. For = if 9+4 (M"), g non- 
negative integer, and A n = y/JI^, we assume that 

^2\ P \\e p \\v < oo. (5) 

pGN 

The process , x) is then formally defined by 

B H (t,x) := y/QY,ep(?)P? (*) = E Vp(^)^ (*)■ 

The sum is normally convergent in C°' 7 (0, T, V), P almost surely for < 7 < H. Indeed, 
in the same fashion as [17], let 

K (u) = ^2 x P \\ e p\\v\\Pp(-,^)\\c^(o,T) 
so that by monotone convergence 

EK = J2K\K\\vn\f3p\\cOn {0 ,T). 

peN 
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According to [21] Lemma 7.4, for every T > and e > 0, there exists a positive random 
variable f] £ ,T,p where E{|77 £i r )P | 9 } is finite for 1 < q < oo and independent of p since the 
f3p are identically distributed, such that \(3p(t) — 0p(s)\ < r] £ T ^\t — s\ H ~ e almost surely. 
Hence, thanks to (jSJ and picking 7 = if — e, we have EfT < 00, 

< 00, P almost surely, (6) 

and B H defines almost surely an element of C°' 7 (0, T, V). As a contrast, a cylindrical 
fractional Brownian motion is defined for a positive self-adjoint Q, which does not 
provide us with almost sure bounds on B H in C°' 7 (0, T, V). Suppose indeed that yfQ is 
a convolution operator of the form ^fQu = g*u for some smooth real- valued kernel q and 
that (e p ) pe N is a real- valued basis of L 2 (R n ). Then, the resulting correlation function is 
stationary (this follows from the convolution and is motivated by ([2])) and 

E{(B H (t, x) - B H (s, x)) 2 } = \t- s\ 2H J> * e p {x)) 2 = \t- s\ 2H \\g\\ 2 L2 

pen 

so that B H belongs to C°- 7 (0, T, L°°(M n , L 2 (fi))) for < 7 < H. As explained in the 
introduction, we are not able to handle such a noise since integration in the probability 
space is required beforehand in order to get some estimates. This is not an issue in the 
context of standard Brownian motions or additive fractional noise, but leads to technical 
difficulties here. 



Fractional stochastic integration. We follow the approach of [T7J [21] based on 
the work of Zahle [32] and introduce the so-called Weyl derivatives defined by, for any 
a G (0,1) and t G (0,T): 



£>oV(0 

D£_f(t) 



1 



r(i -a) 

(-1)° 



/(*) 



a 



fit) - f(s) 



ds 



r(i - a) v(^-^) Q 



/(*) 



a 



(t - s)^ 1 
T f(t)~f(s 



(s-t) 



0+1 



-ds 



whenever these quantities are finite. Above, T stands for the Euler function. Following 
[32], the generalized Stieljes integral of a function / G C°' A (0,T) against a function 
g G C°' 7 (0, T) with A + 7 > 1, A>a and 7 > 1 — a is defined by 



/ fdg := (-l) a / D^ + f(s)D 1 T Z a g T 4s)ds 1 
Jo Jo 

with gT-(s) — g(s) — g(T—). The definition does not depend on a and 



(7) 



fdg :-- 



fHo,t)dg. 



The integral can be extended to different classes of functions since, see 



fdg 



< 



a,iA a (g), 
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, \g(t)-g(s)\ , f W)-g{s)\ A 

A a (g) := ^ sup — — — + a / — — dr 



where 



r(l - a)T(a) <s<t<T V (t-s) 1 -* J s (r - s) 2 ~ a 

so that the integral is well-defined if / G W a> i(0, T) and A a (g) < oo. Besides, the 
fractional integral satisfies the following change of variables formula, see [32]: let F G 
C\R x [0,T]), g G C°' A (0,T) and d x F{g{-),-) G C°' 7 (0,T) with A + 7 > 1, then 

F(g(t),t)-F(g(s),s)= [ d 2 F(g(r),r)dT + f d 1 F(g(r),T)dg(r), (9) 

where dj-F, j = 1, 2 denotes the partial derivative of F with respect to the j coordinate. 

For some Banach space U and an operator-valued random function F G W atX (0, T, C(V, U)) 
almost surely for some a G (1 — H, |), the stochastic integral of F with respect to B H 
is then formally defined by 



C F s dB? :=J2 X p f Fs(e p )d^(s) 
P eN Jo 



(10) 



The integral defines almost surely an element of U for all t G [0, T] since by Jensen's 
inequality for the second line 

/ F s {e p )d^{s) < ^A p A a (/3f)||||F s ( 

e p) \\a,l\\u 

pGN ^° u peN 

< C||-Fs||a,l,£([/,V) X p\\ e p\\vAa((3p) 
p€N 

and as shown in |17j . 

A p ||e p || vA a ({3p(-, oj)) < 00 P almost surely. (11) 

peN 

Hence (TTQ]) is well-defined and the convergence of the sum has to be understood as the 
P almost sure convergence in U. 

We will use the following two results: the first Lemma is a generalization of the 
change of variables formula (Q to the infinite dimensional setting, and the second a 
version a the Fubini theorem adapted to the stochastic integral. Their proofs are given 
in the appendix. Below, V = H q+i (R n ). 

Lemma 2.1 Let F : V x [0, T] — > C be a continuously differentiable function. Let d\F 
be the differential of F with respect to the first argument and d 2 F be its partial derivative 
with respect to the second. For every v G V and B G C°' 7 (0,T, V) for any < 7 < H, 
let (j){t) := diF(B t ,t)(v). Assume that G C°' A (0,T) with A + 7 > 1, and that there 
exists a constant Cm > such that, for all B with \\B\\ c o,-i(o,t,v) < M: 



le°. A (o,T) < CatIMIv- (12) 
Then, we have the change of variables formula, V(s, t) G [0,T] 2 , P almost surely: 

F(B*, t) - F(Bf, s) = f d 2 F(B?, r)dr + J] A p f d 1 F(B?, r) (e p )d^(r). 

Js p( - N Js 
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Lemma 2.2 Let F G W a ,i(0,T, C(V, L^IR™))) with 1 - H < a < \. Then we have: 



„ cN Js \J]R™ J JW 



F T , x dB^ ) dx. 



3 Main result 

We present in this section the main result of the paper. We precise first in which sense 
(HJ is understood. We say that G C°(0, T, H q (R n )) n C°' 7 (0, T, if«~ 2 (R n )), for all 
< 7 < H , g non-negative integer, is a solution to ([!]) if it verifies for all test function 
w G C\0, T, H q+2 (R n )), for all t G [0, T] and P almost surely 

mt),w(t)) ~ (*o, w(0)) = f (s), 9 S «;( S )) rf S 

Jo 

-i/ (y(s),Aw(s))ds + i [ (V(s),w(s)dB?) +i [ (g(V(s)),w(s)) ds, (13) 
Jo Jo Jo 

where the term involving the stochastic integral is understood as 

/' 

Jo 



The latter is well-defined since the mapping F s : e p i-> (We p , it?) belongs to G C°' 7 (0, T, C(V, R)) 
thanks to standard Sobolev embeddings for n < 3. We assume the following hypotheses 
on the non-linear term g: 

H: We have g{e lB ^ x ^) = e^ t,x > g{^) for all real function 9, and for any ^1,^2 m 
H q (M n ) with H^illifa < M, % — 1,2, there exist p G {0, • • • , q} and positive constants 
Cm and such that 

||^(*i) -g(y 2 )\\HP < C' m \\^ 1 -^ 2 \\hp- 

The main result of this paper is the following: 

Theorem 1 Assume that H is satisfied. Suppose moreover that (j3j) is verified for 
V = if 9+4 (R n ), q non-negative integer. Then, for every \I/o £ H q (W 1 ), there exists 
a maximal existence time Tm > and a unique function ^ G C°(0,Tm, H q (W n )) n 
C^^Tm,/^- 2 ^)), < 7 < H, verifying (JT3J) for all t G [0,T M ] P a/mosi surely. 
Moreover, admits the following representation formula: 

q?(t) = e~ lB "u(t, 0)^o + e" ii3f f U(t, s)e lB " g(V(s))ds, (14) 

Jo 

where U = {U(t,s)} is the evolution operator generated by the operator 

iA B n = ie iB " o A o e~ il 



If in addition = 0, then for all t G [0, Tm] the charge conservation holds: 

||*(t)|| L2 = ||*(0)|| L2 . 
If g is globally Lipschitz in H q (M. n ), then the solution exists for all time T < oo. 

When d = 3, a classical example of a non-linearity satisfying H for q = p = 1 is 
g{ty) = V[ty]ty, where V[*f>] is the Poisson potential defined by 

Jr3 \x-y\ 

Indeed, g is locally Lipschitz in if^M 3 ): let \I/i,\I/ 2 G i^QR™); thanks to the Hardy- 
Littlewood-Sobolev inequality [23], Chapter IX. 4, as well as standard Sobolev embed- 
dings, we have 

\\w[*i] - vy 2 [* 2 ]|| L 3 < cih^i 2 - |^ 2 | 2 || L 3 < cw^ - * 2 \\v>\\Vi + * 2 |ki 

and direct computations yield 

< C||*i|li2 + C|pi| 2 || L2 . 

Hence, 

\\g(*i) - g(V 2 )\\m 

< C\\V[^ 1 ]\\ L oo\\^ 1 - y 2 \\ L2 + C\\^ 1 - y 2 \\ H i\\^ 1 + ^Ww-W^zWv 

+ C||^[^i]|| L oo||V^i - V* 2 || L 2 + C||*i - ¥ 2 |MI*i + * 2 ||hi||*2||l« 

<c(II*i||hi + II*2||hi)||*i-*2||hi. 

Another example is given by power non-linearities of the form gfi?) = /i|\I/| 2cr \l/ for 
some /j, G R and a > 0. A L°° bound is needed on \l/ for H to be verified. When n > 1, 
we set then q = 2 and obtain, for all a > |: 

\\gm\m < C\\V\\% +1 , Wg^t) - g(*2)\\v < C||* 2 + *i||&||*i - tfzllw, 
while it can be easily shown that H is verified for n — 1 and q — 1 for all a > 0. 

Remark 2 In order to both lower the spatial regularity assumptions on B H , ty Q , g and 
to obtain global-in-time results, it is natural to consider the energy conservation identity 
(derived formally by multiplying (j25J) by dt^p and integrating, and can be justified for 
classical solutions when q > 2 using the regularity of <p of Theorem and Lemma \2.1\) 
that reads for g = for simplicity: 

hv*(t)\\v = Jl|V*o||| 2 - 9f / / !(i)W(s) ■ VdB^dx. 
£ * Jo Jm. n 

Unfortunately, it is not clear to us how this identity can be used in order to obtain 
estimates on \\ V^\\w a 1 (o,t,l 2 ) f or 1 — H < a < ^ that would depend only on ||V\I/o||l 2 
and \\B H \\ c on(o,T,w l 00 )! \ < 7 < H. Indeed, following the lines of the stochastic ODE 
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case of I21\l in order to treat the stochastic integral and use the Gronwall Lemma, what 
can be deduced from the above relation is an estimate of the form 



\Vnr)\\h\\ WaM o, T) <C + C [ f(s)\\V*(s,.)\\* ds 



for some positive integrable function f and where the constant C depends on ||V\I/o||l 2 
and \\B \\co>i(o,t,w 1o °)- This does not yield the desired bound since we cannot control 
the term II V\I/(s, -)ll 2 „ by II II V\l/(s, OII^IL, . Hence, as opposed to the stan- 

W al {0,T,L 2 ) M" \ '"L \\W a ,l(0,l) 

dard Brownian case, energy methods do not provide us here with an H 1 global-in-time 
estimate. 

Remark 3 When q > 2, then ^ is a classical solution to ([!]) in the sense that it satisfies 
for allt G [0,T m ], P a.s., x a.e.: 

tf(i) = tf(0)+i / M>{s)ds-i [ V(s)dB B -i [ g(^{s))ds. 
Jo Jo Jo 

A proof of this result is given in the appendix. 

The rest of the paper is devoted to the proof of Theorem [TJ The starting point is to 
define <p(t,x) = e iB " {t > x) m(t, x), to use the invariance of g with respect to a change of 
phase and to formally apply Lemma 12.11 to arrive at 

id t ip= -A B Hip + g(ip). (15) 

Remark that A b h can formally be recast as 

A b h = A - 2iVB? ■ V - |V£f | 2 - iAB?. 

In section HI we construct the evolution operator U = {U(t, s)} generated by %A b h and 
obtain the existence of a unique solution to the latter magnetic Schrodinger equation. 
In section [5j we use the regularity properties of the function <p together with Lemma 
12.11 to prove that \I/ = e~ tBt tp is the unique solution to ([T]). The existence follows 
from showing that e~ lBt tp is a solution to ( fTBT) . The uniqueness stems from a reverse 
argument: owing a solution \1/ to ( fl3l) with the corresponding regularity, we show that 
tye lBf is a solution to (fl5l) . This requires some regularization since the function e~ lBt z 
for z smooth cannot be used as a test function in (fT3|) . as well as the interpretation of 
a classical integral involving a full derivative as a fractional integral. 



4 Existence theory of the magnetic Schrodinger equa- 
tion 

The first part of this section consists in constructing the evolution operator U. We 
follow the classical methods of Kato (16] and [22]. The second part is devoted to the 
existence theory for the linear magnetic Schrodinger equation, which is then used for 
the non-linear case. 
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4.1 Construction of the evolution operator 

We follow here the construction of [22], Chapter 5. Let X and F be Banach spaces 
with norms || • || and || • ||y, where F is densely and continuously embedded in X. For 
t G [0, T], let A(t) be the infinitesimal generator of a Cq semigroup on X. Consider the 
following hypotheses: 

(HI) {y4(t)}f € [ 0i T] is such that there are constants u and M > 1, where }u , oo[c p(A(t)) 
for t G [0, T], p(v4(t)) denoting the resolvent set of A(t), and 



it 



< Me wo E ^=! % s,- > 0, < *! < t 2 < • • ■ < T. 



(H2) There is a family {Q{t)}te[o,T\ of isomorphisms of F onto X such that for every 
y G F, Q(t)v is continously differentiable in X on [0,T] and 

where C(t), < t < T, is a strongly continuous family of bounded operators on 
X. 

(H3) For t e [0,T], F C D(A(t)), is a bounded operator from F into X and 
t is continuous in the C(Y,X) norm. 

We then have the following result, see [T6], or p2j, Chapter 5, Theorems 2.2 and 4.6: 

Theorem 4 Assume that (H-l)-(H-2)-(H-3) are verified. Then, there exists a unique 
evolution operator U = {U(t, s)}, defined on the triangle Ay : T > t > s > such that 

(a) U is strongly continuous on A? to C(X), with U(s,s) = I, 

(b) U(t,r)U(r,s) = U(t,s), 

(c) U(t,s)Y C Y, and U is strongly continuous on At to C(Y), 

(d) dU(t,s)/dt = —A(t)U(t, s) , dU(t, s)/ds = U(t, s)A(s) , which exist in the strong 
sense in C(Y,X), and are strongly continuous At to C(Y,X). 

In the next result, we show that for suitable functions B, the operator = ie lB o A o 
e -iB g enera tes an evolution operator U. 

Proposition 5 LetX = L 2 (WL n ) andY = H 2k (R n ), k>\, and MB e C°(0, T, H 2h+2 (R n )). 
Then, the operator iAg generates an evolution operator U satisfying Theorem^ and U 
is an isometry on L 2 (M. n ). 

Proof. We verify hypotheses (H-l)-(H-2)-(H-3) for A(t) = iA B . Let A Bt := A + L(t) 
with 

L(t) = -2iVB t ■ V - \VB t \ 2 - iAB t . (16) 

First, for t fixed in [0, T], the Kato-Rellich theorem [24] yields that A Bt is self-adjoint on 
D(A) = H 2 (W). Indeed, using the regularity B e C°([0, T], H 4 (U. n )), it is straightfor- 
ward to verify that L(t) is symmetric and A-bounded with relative bound strictly less 
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than one. We also obtain that D(A Bt ) = H 2 (R n ), Vt G [0, T}. Stone's theorem [23] then 
implies that for t fixed, iA Bt is the generator of a Co unitary group on X. Moreover, 
— A Bt is positive, so that the spectrum of iA Bt lies in i[0,oo). We therefore conclude 
that the family {jAbJ^^t] satisfies hypothesis (H-l). 

Regarding (H-2), let Q = A^ + J, where / is the identity operator and 

n 

A(*) = (-i)*Z)^' fc ^ L 

3=1 

The operator Q is a positive definite self-adjoint operator on H 2k (R n ), and an isomor- 
phism from H 2k (R n ) to L 2 (IR 2 ). It is also obviously continuously differentiable since it 
does not depend on t. Moreover, 

QAstQ- 1 =A Bt + iQ,A Bt ]Q-\ 

where [A, B] denotes the commutator between two operators A and B. We have the 
following Lemma: 

Lemma 4.1 Fork > 1, let B G C°([0, T], H 2k+2 (R n )). Then [Q, A^Q- 1 G £(L 2 (M n )). 
Proof. We have [Q, A Bt ]Cj _1 = [A( fc ), and using the product rule 

[A( fc) ,L(i)] = 

E E ( ^ ) (-^{VdJlW • - {^|V5| 2 }^ - i{AE$ZB}8r,) 

j=l p=0 N ' x ' 

where ( 2 p fc ) is the binomial coefficient and there are as usual no terms corresponding to 
p = 2k because of the commutator. Using Standard Sobolev embeddings for H 2k+2 (M. n ) 
when n < 3, we have for j — 1, . . . , n that d 2k t+xB G LflP x for p = 6, p < oo and p = oo 

when n = 3, 2, 1, respectively, and d^ q B G L^L^ for g < 2fc. Together with the fact 

that Cj -1 is an isomorphism from L 2 (IR n ) to if 2fc (IR n ) C W 2k ~ 2 '°°(R n ), this is enough to 
insure that [Q, A^Q" 1 G £(L 2 (M n )). □ 

Hypothesis (H-2) is then verified with C(t) = [Q, A Bt \Q _1 , the strong continuity of 
C following from the continuity of B. 

Finally, (H-3) follows easily from H 2k (R n ) C D(A Bt ) = H 2 (R n ), k > 1, and that 
.B G C°(0,T, H 2k+2 (R n )). We can thus apply Theorem H] and obtain the existence of an 
evolution group U generated by iA Bt . The fact that U is an isometry on L 2 (R n ) is a 
consequence of 3t,i(A Bt {p, <p) = for every ip G if 2 (IR n ). □ 

Remark 6 When B G C 1 (0, T, i7 2 (IR n )), a classical choice JM^forQ is Q(t) = XI-A(t) 
for A in the resolvent set of A. This allows to lower the spatial regularity of B but is 
not verified when B = Bf . Notice that in the case when B = Bf, Proposition^ can 
likely be improved in terms of the required spatial regularity of B since we have not used 
the Holder regularity in time of B^ at all. 
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4.2 Application to the magnetic Schrodinger equation 

We apply now the result of the preceeding section to the differential equation 

d t u = iA B u + /, < t < T, u(0) = v, (17) 

where A B = e iB oAoe~ lB . As for CQ), we say that u G C°(0, T, i/ 9 (lR n ))nC 1 (0, T, i^ 9 ~ 2 (M n )), 
g non-negative integer, is a solution to ( fl7|) if it verifies for all w G C 1 (0, T, if 9+2 (M n )), 
for all i G [0,T]: 

(u(t),«;(t))-(T;,iw(0))= / (u(s), 



(u(s), A B w(s)) ds + i / (f,w(s))ds. 



We have the following result: 

Proposition 7 Lei £? G C°(0, T, i? 9+4 (R n )) ; g non-negative integer, and denote by 
U the evolution operator of Proposition Then, for every v G H q (W l ) and f G 
C°(0, T, #"(M n )) ; i/je function 

u(t) = U(t, 0)v + f U(t,s)f(s)ds (19) 
./o 

&e/ongs to C°(0, T, Lf 9 (M n ))nC 1 (0, T, iJ 9_2 (IR n )) and f/ie wnigwe solution to (fEZ} . More- 
over, it satisfies the estimate, for allt G [0,T]: 

IKi)b«<C|M| ff? + C f \\f(s)\\ Hq d Sl (20) 

Jo 

where the constant C depends on ||-B||co(o ) T,fl'9+ 4 (K n )) w^en g ^ 0. 

Proof. Consider first the case q = 2k with fc > 1. The result then follows from 
|16j . Theorem II and the equation (jT7|) in order to obtain the regularity on d t u. The 
cases q = 2k — 1, > 1, and g = are treated by approximation: choose for instance 
sequences £? £ G C°(0, T, Lf 9+9 (IR n )), w e G # 9+5 (M n ) and / e G C°(0, T, H q+5 (R n )) such 
that as £ —7- 0: 

£? £ —>£? in C (0,T,if 9+4 (R n )) (21) 
w £ -»■ v in H q (R n ) (22) 
/ £ / in C°(0,T,#V))- (23) 

Applying the result when q = 2k with fc > 1, the corresponding smooth solution u £ to 
(P3J when g = 2fc - 1 belongs to C°(0, T, # 2fc+4 (M n )) with <9 t w e G C x (0, T, i? 2fc+2 (IR n )), 
with the convention that k = | when g = 0. In order to pass to the limit, it is proven 
in [16], Theorem V, that if iA Be converges to iA B in C(H 2 (M. n ), L 2 (W n )) a.e. t, and 
\\A Be \\c(H 2 (R n ),L 2 (R n )) is uniformly bounded in t independently of e, then 

U e (t,s)->U(t,s) in C(L 2 (R n )) uniformly in (t, s), (24) 
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where U is the evolution operator associated to B. These latter conditions are direcly 
satisfied because of (12T|) . We then write: 

u £ (t) = U e (t,0)v E + [ U £ (t,s)fe(s)ds, Vte[0,T] 

J 

= U(t,0)v+ [ U( y t,s)f( y s)ds + R 1 £ +R 2 £ =u + R 1 £ +R 2 £ 
Jo 

Rl = U £ (t,0)(v £ -v)+ [ U £ (t,s)(f £ (s)-f(s)) 

Jo 

R 2 £ = (U £ (t,0)-U(t,0))v + [ (U £ (t,s)-U(t,s))f(s))ds. 

Jo 

Using ( |24|) and the strong convergence of v £ and f £ , we then obtain that u £ — > u in 
C°(0, T, L 2 (R n )). Assume first that q 7^ 0. In order to get the announced better regu- 
larity on u, we use the fact that u £ G C°(0, T, C 2k+2 (M. n )) and <%w £ G C^O, T, C 2k (R n )) 
thanks to standard Sobolev embeddings for n < 3. We can then differentiate equation 
(JSJ), and find using the representation formula 

D^u £ (t) = U £ (t,0)D^v £ + [ U £ (t,s)(D^f £ (s) + [D f} ,L £ (s)}u £ (s))ds, (25) 

Jo 

where 1 < \j3\ < q and 

D ^' = d^ x '" x &p ^ = " ' ' 1/31 = A + ' ' ' + ^' 

and L e (s) is defined in ( fT6l) with i? replaced by B £ . Only the term involving the 
commutator requires some attention. Using (12ip . we can show that for all s G [0,T], 

||[^,L e ( S )K(s)|| L2 <C|k( S )|U,, 

where the constant C is independent of e. Together with (|2T]l - p2"]l - (p3]l -p4])-(p5]) and 
the Gronwall lemma, this yields a uniform bound for u £ in C°(0, T, H q (M. n )). Using 
this latter bound along with (I21]l - (j22]l - (j25)l - (j21]> - (125)1 and equation (H7J) for the smooth 
solution u £ in order to estimate dtu e , it is then not difficult to show that (u e ) e is a Cauchy 
sequence in C°(0, T, H q (R n )) n C 1 (0, T, if 9 ~ 2 (IR n )), whose limit u satisfies estimate flU 
and (jTSI) . When q — 0, it suffices to use equation (JT7J) for the smooth solution w e in 
order to show that {d t u £ ) £ is Cauchy in C°(0, T, if ~ 2 (IR n )). This proves the existence, 
the representation formula ( 1T91) and estimate ( 1201) . 

Uniqueness is straightforward in the case q > 1 since solutions to (1171) are regular 
enough to be used as test functions and to obtain after an integration by part that 
3(V(e~ lSt ii), V(e _lB 'u)) = 0. When g = 0, we use the adjoint formulation of (TIT]) . 
The difference between two solutions to ([171) satisfies in the case of a test function 
w G C°(0,T, ff 2 (M n )) nC^T, L 2 (M™)), 

(u(t),u>(t)) = / (it(s),9 s w(s) + (zA Ss )*w(s))ds, 
•/ 
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where (iA Bs )* = -iA Bs is the adjoint of iA Bs . Let t G [0,T], pick some w G if 2 (IR n ) 
and let w(s) = z(t — s) where z(s) is the solution to d s z(s) = (iAB t _ s )*z(s), z(0) = wq, 
< s < t. Adapting Proposition to the operator {iA Ba )* , Theorem 0] yields that z G 
C°(0, T, # 2 (M n )) n C^O, T, L 2 (R™)). Hence, <9 s w(s) + (zA B J*w(s) = z a.e., w(f) = w 
and it comes, for all t G [0, T\. 

(«(t),w o ) = 0, V«i Gff 2 (l n ), 

so that it = 0. This ends the proof. □ 

We use the result of the last Proposition to prove that the non-linear magnetic 
Schrodinger equation 

dt<p = iA B H<p + g(<p), 0<t<T, u(0) = tt o , (26) 

admits a unique solution P almost surely in the same sense as ffl8|) : 

Theorem 8 Assume that H is satisfied. Suppose moreover that (JSJ) is verified for 
V = H q+A (W n ), q non-negative integer. Then, for every ty G // 9 (lR n ), there exists 
a maximal existence time T M > and a unique function <p G C°(0, T M , H q {R n )) R 
C 1 (0,T M ,H q - 2 (R n )) verifying ([26]) for t G [0, T A/ ] w/i«c/i arfmrfs the following represen- 
tation formula: 

<p{t) = U(t,0)V + I U{t,s)g{cp{s))ds, 
Jo 

where U = {U(t,s)} is the evolution operator generated by the operator 

iA B n = ie iB " oAo e~ iB " . 
If moreover Qg((p)Tp = 0, then for all t G [0, Tm] 

h(t)\\L* = \\m\\v. (27) 

If g is globally Lipschitz on H q (M. n ), then the solution exists for all time T < oo. 

Proof. The proof is very classical and relies on Proposition [7] and a standard fixed 
point procedure. First of all, (JSJ) insures that P almost surely, B H G C°(0, T, H q+A (W 1 )), 
which allows us to define an evolution operator U according to Proposition The rest of 
the proof follows the usual arguments of for instance [22], Theorem 1.4, Chapter 6, that 
we sketch here for completeness. Given \l/o i n H q (M. n ) and (p G C°(0, t\, H q (M. n )) for some 
t x > to be fixed later on, denote by u := F(ip) the solution to (HZj) in C°(0, t u H q (R n ))n 
C 1 (0,fi,if 9 - 2 (M n )) where / = g(<p) belongs to C°(0, t x , H q (R n )) thanks to hypothesis 
H. Using the latter, estimate f!2Up . following the aforementioned theorem of [22], one 
can establish the existence of M > and a time t\(M) such that F maps the ball of 
radius M of C°(0, ti, H q (M. n )) centered at into itself. In this ball, the function g being 
uniformly Lipschitz on H p (M. n ), < p < q according to hypothesis H, existence and 
uniqueness of a fixed point of F in C°(0, t\, H q (W 1 )), denoted by tp*, follows from the 
contraction principle. Moreover, this solution verifies the representation formula f[T9l) 
with / = g(tp*) G C°(0, tx, H q (M. n )) and according to Proposition [TJ belongs in addition 
to C 1 (0, ti, H q+2 {R n )) and satisfies f!2B|) . The existence of a maximal time of existence 
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Tm is established following the same lines as [22]. When g is globally Lipschitz in 
H q (M. n ), then Tm < oo by the Gronwall Lemma. 

Regarding the conservation of charge (1271) . the case q > 1 is direct since the solution 
ip* is regular enough to be used as a test function in ( TT8"j) (after interpretation of (•, •) 
as the H^ 1 — H l duality pairing when q = 1) and it then suffices to take the imaginary 
part of the equation. When q = 0, we use a regularization procedure very similar to 
that of the proof of Proposition [7J the details are left to the reader. □ 

5 Back to the stochastic Schrodinger equation 

We apply the result of the last section to prove Theorem [TJ Owing the solution (p of 
Theorem [HI it suffices to show (i) that e~ lBt ip is a solution to f [T3|) . which will follow 
from the regularity of <p and Lemma I2TTI this yields existence; and (ii) that all solutions 
to f[T3"j) with the corresponding regularity read e~ lBt u where u is a solution to fT2"6"|) . 
which yields uniqueness since ([26]) has a unique solution. As explained earlier, the last 
step requires a regularization procedure since test functions of the form w = e~ lBt z, 
with z smooth, are not differentiable in time and cannot be used directly in (fT3l) . In 
the whole proof, T denotes some time T > T M . 

Proof of Theorem [TJ. Existence. Let <p be the unique solution to f[26|) according to 
Theorem [8] and define, for any test function w G C 1 (0, T M , H q+2 (M n )), 

F(B?,t) = (e-* B ?<p(t)Mt))- 

We verify that F satisfies the hypotheses of Lemma 12.11 First of all, F is clearly 
continuously differentiable w.r.t. the first variable and for all v G H q+4 (M. n ), let (p(t) : = 
d 1 F(B^ I ,t)(v) = i(e~ iB ?(p(t)v,w(t)). Second of all, we need to show that <p G C°' A (0,T) 
for A verifying A + 7 > 1, together with the bound (TT2"1) . To this goal, we have for 
(t,s) G [0,T M ] 2 : 

<Kt)-4>(a) = i (V^ -e^')<p(t)vMt)) +* (e^W) -<p(s))vMtj) 
+i (e" iB "(p(s)v, (w(t) - w(s))\ 
:= T 1 + T 2 + T 3 . 

We treat each term separately. We have, using standard Sobolev embeddings for n < 3: 
\Ti\ < C\\<p(t)\\ L 2\\w(t)\\ L o a \\v\\ L o \\Bf -B?\\v 

< C(t - sy\\v\\ H q+*\\B? ||c0.7(0,T M ,L2) 

< C(t-sy\\v\\ Hq+ 4, (28) 

for all < 7 < H. Regarding the term T 2 , notice that the product e tB "vw belongs to 
H q+2 (R n ) when n < 3, so that since d t cp G C°(0, T M , H q ~ 2 (R n )), we can write 

(e- lB "(tp{t) -<p{s))v,w(tj) = f (d TV (T),e tB "vw(t)) Hq -2, Hq+ 2dT, 
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where when q > 2, the pairing (•, -) Hq -2 Hq +2 is replaced by the L 2 inner product. Hence, 
\T 2 \ <\t- s\\\w(t)\\ Hq +2\\v\\ Hq +4Bf \\c°(o,T M ,m+*)\\dM\eo(o,T M ,m-^ < C\t- s\\\v\\ Hq +4. 

Estimation of T 3 is straightforward and leads to a similar estimate as above. This, 
together with f[2"gj) yields that 

ll0llc°^(o,r M ) < C|MIh9+4. 

Since \ < H, we can pick H = | +e and 7 = | + 1 such that 27 > 1 and the assumption 
on of Lemma 12.11 is verified. It remains to show that 82 F exists and is continuous, 
and this is a consequence of the fact that d t (f G C°(0, Tm, H q ~ 2 {R n )). Applying Lemma 
12.11 then yields 

(e- w "ip(t), w(t)) - w(0)) = f\dMr), e w "w(r)) Hq -2 Mq+2 dT (29) 

Jo 




In order to conclude, picking w(t,x) = w(x) G H q+2 (M. n ) in (Tl8j) with / = g(<p), it 
comes that (126|) is verified in H q ~ 2 (M. n ) for all t G [0, Tm] and almost surely. This yields 
d T (p = iA B H(p — ig(ip) in H q ~ 2 (M. n ), and replacing d T (f by its latter expression in (129|) . 

setting ^ = e - iB **V G C°(0,T M , H q (R n )) n C°- 7 (0, T A/ , i/^ 2 (M n )) for all < 7 < 
finally yields ([13]). 

Uniqueness, Step 1: regularization. Starting from a solution ^ to ( fl3l) with the 
above regularity, we would like to choose the test function w = e~ lBr z for some regular 
function z in order to recover the weak formulation of (1261) . which admits a unique 
solution. This is not allowed of course since B H is not differentiable. The solution is to 
use the Holder regularity of \1/ in order to reinterpret the term 

/ {^{s),d s w{s))ds 
Jo 

clS 8b fractional integral. To this end, let 

B H >%t,x) :=^A p e p (x)/3f e (t) 

pen 

where (3p' e is a C l regularization of (3^ such that (3p' e — > (3p in C 0,7 (0,T) almost surely 
for all p and H/^jf^Hc^to.T) < ||Af \\c°<f{o,T), < 7 < H. We have 

B H,e^ B H in c°' 7 (0,T,if 9+4 (R")), P almost surely. (30) 

Indeed: 

\\B H ' £ — B H \\ c o,i(o t T,Hi+ 4 ) < ^ ^pll e pl| J H""+ 4 l|/5p I/ ' e — /^Hc.T^T) 

and 

||/3^' e - /3jf ||c°.7(o,t) < 2||^||cO,t(o,T), 
17 



which, together with the convergence of (3p' e to j3p in C ' 7 , ([6]) and the Weierstrass rule 
gives the desired result. Set w £ = e~ iB "" z where z E C 1 (0, T, H q+2 (R n )). Then 

/ (y(s),d s w E (s))ds 
Jo 

y(s)e w "'\d s z(s)^ ds-iJ2K I (^{s)e iB "'\ z{s)e p ) (^ £ (s))'ds, 

where all permutations of sum and integrals were permitted since the series defining 
B H,£ is normally convergent in C 1 (0, T, H q+A (W n )). We now use the fact that for a 
continuously differential) le function /, Df_f — > f as a — > 1, see [32] section 1, and that 
D^ 13 f = D^_D^_f, where the operator D^_ is defined in section EJ We then introduce, 
for \ < 1 — [i < 7, 

I £ P ■= f (^(s)e^\e p ) (^(s) - 0?*(t))'d8 

= lim I £ p a := lim / (*(s)e iB ?' S , e^( a )) A-"^t(/?f' £ )*-(«)^> 

where (/3f- £ ) t _(s) = - /3f' e (r). Owing the fact that (*(s)e iB "' E , z(s)e p ) G 

C°' 7 (0, Tjvf) for any < 7 < H and using the fractional integration by part formula of 
[32] section 1, we find 

IT = (-I)*"" jTzftr" (^(s)e^,z( S )e p ) £ ) t _( S )d S . 

Moreover, we can send a to one above thanks to dominated convergence since the term 
(ty(s)e lBa ' , z(s)e p ) belongs to C 0,7 (0, T M ) and 1 — [i < 7 in order to obtain 

We derive below some estimates needed to pass to the limit £ — > 0. 

Uniqueness, Step 2: uniform estimates. Recall that w £ = e~ lBt ' z and define first, with 

w = e * z: 

(j) £ (s) = (^(s),(w £ (s)-w(s))e p ) := (V(s),r £ (s)e p ). 
In order to estimate JJ o 7 M e , we write 

fit) - 4>%s) = (tt(t) - V(s),r £ (t)e p ) + (r e (t) - r e (s))e p ) := 7\ + T 2 . 

For the term T 1; we use the C 0,7 regularity of \& in H q ~ 2 , while we use that of r e for T 2 . 
It comes with the help of standard Sobolev embeddings: 

1^1 1 < \t - s\ J \\^>\\cO,i(o,T M ,Hi- 2 )\\l r £ {t)\\Hi+ 2 \\e p \\Hi+ 2 

\T 2 \ < \t - s| 7 ||^(t)|| L 2||ep|| L °o||r e ||co.7(o,r,i 2 )- 
Since 1 — fi < 7, this gives 

11^11^1-^,1(0^/) — C\\e p \\ H g+4\\B H,£ — B H \\ C 0,j(Q jT>H q+4Q R n)y (31) 
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On the other hand, using the notation of section [21 we find 

\Dl(^) t _(s)\ < Ai_ M (/3f £ ) < C\\^\\ c o,- i{0 ,t) < C\\ff\\o ivn < oo, (32) 

and 

|£>t(/f' £ - #V(*)I < Ai- M (/?f £ - ) < C||^ e - /3f || c o.,(o,T). (33) 
Uniqueness, Step 3: passing to the limit. We have all needed now to pass to the limit 
in the weak formulation (EE)). Plugging w £ (t) = e~ iB "' E z E C 1 (0, T, iJ 9+2 (IR n )) yields 

(*(t),u; e (t)) - (*o,*(0)) = [ ms),d s w%s))ds 

Jo 

(*(s),Aw £ (s))ds + i / (^{s),w £ (s)dB B ) +i [ (g(^(s)),w £ (s))ds. (34) 

Jo JO 



We have 



/ (^(s),<9 s w; e (s))cfe = / (^(s)e^' £ ,^z(s)) - z V I 
Jo Jo v / peN , J 



Using (]HT) - (T30T) - fT3TT) - fT32l ) - (T33T) as well as (jSJ), we can pass to the limit in the latter equation 
and obtain that, Vi G [0,7^/]: 

limjf {^{s),d s w £ {s))ds = (^{s)e iB ?,d s z{s)^ds-i J (^{s)e iB ° ' ,z{s)dB^. 

Similar arguments can be employed to pass to the limit in the remaining terms of ([3 
The stochastic integrals simplify and we are left with 

(v(t)e iB ?,z(t)) - (*o,*(0)) = f (^{s)e iB ? ,d s z{s)) ds 



- % J* A(e~ iB " 'z(s)fj ds + i (g{*{s% e~ iBf z{s)) ds. (35) 

Hence, tye tB " verifies the magnetic Schrodinger equation (ITSl) with / = g(^(s))e lB " = 
g(^/(s)e iB "). Since the latter admits a unique solution cp e C°(0, T M , F 9 (M ri ))nC 1 (0, T M , H"~ 2 (R n )) 
according to Theorem El we can conclude that (CQ) admits a unique solution. The repre- 
sentation formula (IT4"|) follows then without difficulty with the identification ^e lBt = tp 
and Theorem El This ends the proof of Theorem [TJ 



6 Appendix 

6.1 Proof of Lemma |2JJ 

First of all, we know by section [2] that B B belongs to E :— C 0,7 (0, T, V), P almost surely 
for < 7 < H. We proceed by approximation in order to apply the change of variables 
formula ([9]) valid in finite dimensions. Let 

N 

B?' N (x) :=£v*(3))flf(*), 

p=0 
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so that, 

B? > N ->■ Bf in £, P almost surely, (36) 

thanks to (J3J) and ([6]). We have moreover the bound \\B?' N \\ E < \\Bf\\ E := M. bmce 
F is C 1 , and e C°' A (0, T) with A + 7 > 1, we can use © and find, for < s < t < T 
fixed: 

F(B?' N ,t)-F(B?> N ,s)= / d 2 F(B?> N ,T)dr + J2KJ d 1 F(B^ N ,r)(e p )d^(r). 

p=0 ^ s 

By continuity of F, it is direct to pass to the limit in the left hand side. The same holds 
for the first term of the right hand side thanks to dominated convergence and the fact 
that d 2 F is continuous and B^' N is bounded in E independently of N . Regarding the 
last term, let <${t) = d x F{B^ N , r)(e p ) and 

fp ■= f ${T)dfi{T) = {-IT f D« +( f>»(T)D}:°<(P*) t „(s)dT, 

by (I7j) for some a verifying a < A and 1 — a < 7. Since H-B^^H^ < M, we have by (fl2l) 
^(r) - $(8))\ < |t - «| A " a " 1 lKllco.*(o,r ) < Cm|t - s^Mv, 



IT — S\ 



where A — a > 0. The latter estimate, dominated convergence, ( 1361) together with 
the continuity of d\F yield first that Df + (j)p (r) — > Dg + <p p (r) a.e. where <p p (r) = 
d x F(B* ,r)(e p ). Then, since |D£r a G0*)t-(s)| < A Q (/3f) < 00, P almost surely, we 
have 

\D^{r)D}z a {^) t -{s)\ <C\r- s^^Wv (37) 

so that dominated convergence implies that f p N ^f P = J* 4 ) p( T ) d/3p(r), for all pel 
Finally, since 

\\f p I < C^p||e p ||y, 



thanks to fl37|) . and moreover (joj) holds, we can apply the Weierstrass rule and conclude 
that P almost surely: 



N 00 „t 



p=0 p=0 

This ends the proof. 
6.2 Proof of Lemma 12.21 

The hypothesis on F show that the integral on the left is well-defined and that 

VA P I ([ F TtX (e p )<to)dP?(T)= lim VA P / f / F r , x (e p )^ (r) := lim I 
Jo \Jr" J N ^°°pTo Jo \Jr" J 



AT- 
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Morever, for 1 - H < a < |, we have that \D]z a {Pp) t -{s)\ < A a (/5f) < oo, P almost 
surely and D% + F t)X (e n ) G L\(0,T) x R n ) since F t;X (e p ) G W tt ,i(0, T, L^R")). Hence, 
using the definition of the stochastic integral and Fubini Theorem, it comes 



in 



(" 1 ) H f> / f r ^ + ^(e P )(r)] [A- a (^)*-(r)] 

it" \ „_ n JO / JR n 



dr I c&e 



Moreover, P almost surely: 



N 



p=0 



\W a ,i(0,T)\\ L i 



< C||^||iy Q ,i(o,r,£(y,Li)) ^ A p ||e p ||yA a (/3p 

so that thanks to (TTTj) . the series defining converges strongly in L l (M. n ) and almost 
surely This yields 

\imI N =[ ( VA P / F ttX (e n )dp*(T) ) dx 



,p=0 



and ends the proof. 



6.3 Proof of Remark [3] 

For q > 2, using the regularity * G C°(0, T A/ , # 9 (M n )) n C°' 7 (0, T M , H q ~ 2 {R n )) and 
picking w G L 2 (M. n ), we can recast ( !T3|) as 

{V(t),w)-(V ,w) = -i f (£9(s),w)ds + i f (^{s),wdB^) + % [ (g(V(s)),w)ds. 

Jo Jo Jo 

Since the mapping F : e p — > ^(s)we p belongs to C°' 7 (0, T M , C(V, we can use 

Lemma 12.21 together with Fubini theorem to arrive at 

{*(t),w)-{V ,w) = 
which yields the desired result. 
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